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BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

The present invention relates generally to the field of radiation detection and 
imaging technology. In particular, the present invention may be used for processing 
exponentially decaying pulses from a scintillation detector or a similar detector for 
radiation particles. 

2. Description of Related Art 

When a radiation particle (gamma ray, neutron, electron, etc.) is detected in a 
scintillation detector, the scintillation detector will emit light, which is then converted 
into an electronic signal by a photosensor (e.g., photomultiplier tube or photodiode). This 
electronic signal can then be received and processed by electronic circuits. In the period 
after radiation hits the scintillation detector, the scintillation light decays exponentially 
with a time constant x (the time when the light level decays to 37% of the onset level), as 
shown in FIG. 1. 
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FIG. 1 shows energy output by two gamma ray particles over time. Since the total 
amount of light emitted by the scintillation detector represents linearly the energy 
deposited by the radiation particle in the detector, the area or integral under the curves in 
FIG. 1 is a measure of the particle energy. As shown in FIG. 1, area 5 and area 10 define 
5 a measure of the particle energy of the gamma ray particles. Furthermore, the initial peak 
in the light level is also proportional to the radiation energy. Hence both the area 5 and 
peak VI in FIG. 1 may be used to measure the energy of the gamma ray or radiation 
particle. Since the area under the curve (integral of light) includes many more light 
signals than the instantaneous peak light level, the integral (the total amount of light 
10 emitted) is generally used to measure the radiation energy. 

As the radiation flux increases, it becomes increasingly likely that the next 
radiation particle may arrive at the detector while all previous events are still emitting 
light (FIG. 2). In this case, the identity of each individual radiation particle will be lost, 
and several particles will merge into one large signal, as shown in FIG. 2. In this case 
15 neither the peak level (VI or V2 of FIG. 1), nor the integral information (area 5 or area 10 
of FIG. 1) can be used to separate or measure the energy of each particle. In these 
situations, the detection system will fail to respond properly because of erroneous 
measurement. 

It is known that it takes a time period of approximately 4% to collect 98% of the 
20 scintillation light from each radiation excitation. Thus, if the next event arrives at time 
t > 4x, the pile-up-energy error on the next event will be less than 2%. Hence, to keep 
pile-up error small, it is desirable to minimize the chance that two events (radiation 
excitations) will occur in a time period less than 4t. Since the time-lapse between two 
events is a random distribution (i.e., the time-lapse between two events is a random 
25 variable) centered about the "average arrival time", it is generally practiced in the prior art 
to operate the detector so that the "average arrival time" is 10 x (4t) = 40x, to lower the 
random chance of having two events coming closer than 4%. With this lOx "head-room", 
the probability that two events will come closer than 4x would be approximately 10% 

-3- 

25075668.1 



(using Poisson statistics), 
shown in Table 1 below: 



The head-room factor as a function of pile-up percentage is 



TABLE 1 

HEAD-ROOM FACTOR AND PILE-UPS 



Head-Room Factor 


5 times 


10 times 


15 times 


20 times 


pile-up / total (%) 


18 % 


10 % 


6.5 % 


5 % 



5 

Thus, a lOx head-room is a reasonable choice, and is generally practiced in the 
prior art. When coupled with a 4x light-collection time (system dead-time), such a prior 
art detector provides a measured-energy error (due to pile-up) of less than approximately 
2% for approximately 90% of the time, and an energy measurement error (energy 
10 resolution) greater than 2% for approximately 10% of the time. This minimum lOx 

head-room (40x) timing requirement means that the maximum detection-rate should be 
less than 1/(40t) for the scintillation detector. 

As the role of gamma cameras expands in positron coincidence imaging, 
radionuclide therapy dosimetry imaging, and cardiac first-pass imaging, there is a need to 

1 5 significantly increase the count-rate capability of gamma cameras because of the large 
photon flux in these procedures. Single-crystal design, conventional detector electronics, 
and traditional Anger-positioning algorithm may all hinder higher count-rate imaging 
because of the pileup of signals of detected events in the detector. In PET cameras, 
especially the type using Anger-positioning detectors (all commercial systems), pileup 

20 may also be a problem. 

The present disclosure permits a scintillation detector system to operate at a much 
higher event-rate (count-rate) by obviating the lOx head-room factor without pile-up. 
The techniques of this disclosure maintain a greater event-rate with little sacrifice in the 
total amount of scintillation light collected, specifically at a 10 times higher radiation flux 
25 with little or no sacrifice in measurement accuracy. If the fraction of scintillation light 
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collected can be reduced (i.e., if a user is willing to compromise measurement accuracy), 
the present disclosure allows detectors to count at count-rates approximately twenty times 
greater than conventional methods. Techniques of the present invention may be termed 
HYPER (High-Yield-Pileup-Event-Recovery). 

5 This disclosure proposes different implementations to execute HYPER 

mathematics to either improve resolution of the total energy (or position) estimation for 
the HYPER method or to provide a simple way (or a faster way) to further increase 
count-rate using a multiple HYPER zone method. This disclosure also shows how to use 
a HYPER circuit in the application of coincidence imaging. 

10 SUMMARY OF THE INVENTION 

The present invention includes an apparatus for signal pile-up prevention, 
comprising a delay circuit for receiving, holding, and passing an incoming signal; a 
computation circuit for determining a weighted value of the incoming signal; a sampling 
circuit for receiving the weighted value. The sampling circuit passes the weighted value 
15 (which may be passed to an A/D converter) upon receipt of a triggering signal, which 
corresponds to receipt of a next incoming signal at a trigger circuit. In an exemplary 
embodiment, the computation circuit may comprise an amplifier, an integrator, and an 
adder. In an exemplary embodiment, the weighted value is a sum of an integrated value 
and an instantaneous value, and may be a substantially constant value. 

20 An apparatus according to the present invention may also include a smoothing 

circuit connected to the circuit adapted to receive the incoming signal The apparatus may 
also comprise a residual subtraction circuit for reducing the weighted value by a residual 
signal value. The sampling circuit discharges said weighted value upon input of the 
triggering signal. 

25 The present invention may be used in connection with nuclear medicine 

applications, such as a PET or gamma camera, and may be used to determine both energy 
and position information. Such an apparatus comprises a plurality of delay circuits, a 
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plurality of computation circuits, and a plurality of sampling circuits, wherein each of the 
delay circuits receives a different incoming signal from a different output of a gamma 
camera. The delay circuit, computation circuit, and sampling circuit comprise a pile-up 
prevention circuit. 

Particular embodiments will comprise a plurality of pile-up prevention circuits, 
and may include a digital signal processor and fast trigger connected to each of the pile- 
up prevention circuits. Such an embodiment may also comprise an inter-zone detection 
circuit connected to the fast trigger and a multi-zone-trigger processor connected to said 
inter-zone detection circuit, capable of centroid averaging. An exemplary embodiment 
may have a plurality of fast triggers. 

A method for preventing signal pile-up according to the present invention may 
comprise: delaying an incoming signal for a preselected time prior to passing the 
incoming signal; computing a weighted value of the incoming signal; and sampling the 
weighted value upon receipt of a triggering signal from the next radiation particle, 
thereby preventing signal pile-up. Computing may include amplifying the incoming 
signal to obtain an amplified signal, integrating the incoming signal to obtain an 
integrated signal, and adding the amplified signal and integrated signal to obtain the 
weighted value. This method thereby creates a variable signal collection time. 

Another method according to the present invention may determine position and 
energy information of incoming signals without pile-up. Such a method may include: 
delaying at least one prenormalized position signal and a total energy signal; computing a 
weighted value for each prenormalized position signal and the total energy signal; 
sampling the weighted value for each prenormalized position signal and the total energy 
signal, upon receipt of a triggering signal from the next radiation particle. In an 
exemplary embodiment, the prenormalized signals and the total energy signal may be 
corrected by subtracting remnant values of all previous signals. 

In yet another aspect, the present invention comprises an apparatus for 
dynamically detecting energy of each one of a plurality of incoming signals received 
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from a detector, without pile-up of previous incoming signals, including: a delay circuit 
connected to receive an incoming signal from the detector and to pass the incoming signal 
from an input to an output of the delay circuit after a time delay; a trigger circuit 
connected to receive the incoming signal from said detector, and for generating a 
5 triggering signal upon receipt of a subsequent incoming signals at the trigger circuit; a 
computation circuit connected to the output of said delay circuit for determining a 
weighted value of the incoming signal; a sampling circuit connected to receive the 
weighted value from said computation circuit, and for passing the weighted value from an 
input to an output of the sampling circuit upon receipt of the triggering signal; and a 
10 residual subtraction circuit connected to the output of said sampling circuit, for 

subtracting a residual signal value corresponding to a residual weighted value of previous 
incoming signals, and for providing an output signal corresponding to the energy of the 
incoming signal. 

Another aspect of the present invention comprises an apparatus connected to a 
1 5 gamma camera for detecting position and energy information of each one of a plurality of 
incoming signals received by the gamma camera, without pile-up of previous incoming 
signals, including: a first delay circuit connected to receive a first incoming signal from 
the gamma camera, and for passing the first incoming signal from an input to an output of 
the first delay circuit after a first time delay; second and third delay circuit arranged like 
20 the first delay circuit to receive, delay, and pass, second and third signals; a trigger circuit 
connected to receive the third incoming signal from said gamma camera, and for 
generating a triggering signal and a timing mark upon receipt of a next third incoming 
signal at the trigger circuit; first, second and third computation circuits, each connected to 
receive an output of a respective one of the first, second, and third delay circuits, and for 
25 determining a respective weighted value for each of the first, second, and third incoming 
signals; first, second and third sampling circuits, each connected to receive a respective 
one of the first, second, and third weighted values, and for circuits passing the respective 
weighted value upon receipt of the triggering signal; and a digital signal processor 
connected to receive the first, second, and third weighted values, and for subtracting 
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residual signal values corresponding to residual weighted values of previous ones of the 
incoming signals, and for providing an output signal corresponding to a position value of 
the first and second incoming signals and an energy value of the third incoming signal 

5 Another aspect of the present invention resides in a method of obtaining energy 

information for each one of a plurality of incoming signals received from a detector, 
without signal pile-up, comprising the steps of: delaying an incoming signal for a 
preselected time; computing a weighted value of the signal after the preselected time; 
sampling the weighted value upon receipt of a subsequent signal; and subtracting a 
10 residual signal value from the weighted value to obtain the energy information. The 

residual signal value may correspond to a residual weighted value of at least one previous 
incoming signal, thereby preventing signal pile-up. 

Yet another aspect of the present invention resides in a method of determining 
position and energy information of a plurality of incoming signals from a detector 

1 5 without pile-up, comprising the steps of: receiving a first and second prenormalized 

position signal and a total energy signal from the detector; delaying the first and second 
prenormalized position signals and total energy signal for a preselected time; computing a 
weighted value for each of the first and second prenormalized position signals and total 
energy signal after the preselected time; and sampling the weighted value for each of the 

20 first and second prenormalized position signals and total energy signal upon receipt of a 
subsequent first and second prenormalized position signals and total energy signal. 

The methods of the present invention may be used to operate gamma-cameras (or 
other radiation detectors) in very high count-rate situations. The present invention 
25 includes the following features: (a) no compromise in measured energy-resolution in low 
count rates; (b) count recoveries and accurate energy measurement even for gamma-rays 
within a pile-up involving multiple gamma-rays; (c) optimal scintillation-light collection 
in very high count-rate situations; and (d) ability to merge with a multi-zone architecture 
to further increase count-rate capability. The present invention includes algorithms that 
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apply to all triggering gamma-rays (it is to be understood that although gamma-rays are 
discussed herein, the present invention applies to all types of radiation detectors), for 
extracting the correct energy and position of every triggering gamma-ray. 

This disclosure also provides ways to improve energy and position resolution for 
5 the HYPER (High-Yield-Pileup-Event-Recovery) method and illustrates how to use that 
method in the application of coincidence imaging. Conventional single-scintillation- 
crystal design coupled with conventional detector electronics and the traditional Anger- 
position algorithm may all hinder higher count-rate imaging because of the pileup of 
gamma ray signals in the detector and electronics. At an interaction rate of 2 million 

10 events per second for Nal(Tl) scintillators, the fraction of nonpileup events is <20% of 
the total incident events. Hence, the recovery of pileup events may significantly increase 
the count-rate capability, increase the yield of imaging photons, and minimize image 
artifacts associated with pileups. The methodology of this disclosure, coined "HYPER," 
solves the pileup problem. One main idea of HYPER may be stated as (a) dynamic 

15 integrating: if no pileup happened, a signal from an event may be integrated for a fixed 
time (e.g., 3 to 4 times the scintillation decay time) to fully collect all the light (signal) to 
maximize detection accuracy (resolution); however, if a pileup is detected, the integration 
of the current signal stops immediately; (b) calculating a weighted value: a weighted 
value is the estimated total energy that includes the energy of the current event and a 

20 residual energy from all the previous events; (c) remnant correction: calculating a pile-up 
free energy from two consecutive weighted values. 

FIG. 20 shows an exemplary diagram of one embodiment of HYPER techniques. 
The energy and pre-normalized position signals from a gamma camera are delayed for a 
fixed period of time. The energy-signal pulse goes to a triggering circuit that generates a 
25 trigger and a timing mark at the leading edge of the event-signal. The computation circuit 
generates a weighted value estimating a total energy or pre-normalized position signal 
that includes the contribution from the current event and the remnant signals from all the 
previous events. The present-event weighted value is sampled and buffered, and two 
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consecutive weighted values are used for the correction of residual signals from all 
previous events. 

One of the more significant discoveries concerning the HYPER method is that no 
matter how many previous events are piling up onto a present event, only the information 
5 from the present event and the last event are needed to subtract the pileup effects of all 
previous events (say, 5-6 events). This two consecutive-event method makes processing 
efficient and fast. The relation of the input signal, weighted value, and the residual signal 
is illustrated in FIG. 21. 

HYPER methodology may be implemented and improved upon in several ways. 
10 Different implementations may use the same procedure and mathematical formulae as the 
basic HYPER method mentioned above in (a), (b) and (c). Different implementations 
have different advantages and trade-offs in signal noise and processing speed. When the 
signal noise is high, a lower noise weighted value for estimating the total energy may be 
extracted in several different ways: 

15 (1) Creating a weighted-sum of an instantaneous value and an integrated value. 

The weighted-sum is the total energy in theory, but it can be noisy because the 
instantaneous value itself has more noise than the integrated value. To reduce the noise, 
the weighted-sum signal may go through an analog low-pass filter, and the filtered 
weighted-sum may be continuously or multiply sampled. Then, a digital filter such as 

20 median filter or a decay filter may be applied to those samples to generated a low noise 
weighted value. 

(2) Weighting the integrated value. To get the total energy (present energy plus 
remnant energy), the integrated value may be dynamically weighted by a factor. 

(3) Dynamically weighting the sum of the continuous-sampled values of the 
25 instantaneous signal. The sum signal of multiple/continuous sampled values of the 

instantaneous signal is the digital integration of the instantaneous signal. This method is 
similar to the method of (2) described immediately above, the difference being digital 
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integration instead of analog integration. The advantage of digital integration is that no 
analog discharge is necessary, which thereby eliminates analog-discharging dead-time of 
the system. 

To further increase count-rate capability, multiple HYPER zone circuits may be 
5 implemented. Methods are presented in this disclosure to identify a zone that is hit by a 
real gamma. FIG. 25 shows a large detector divided into four zones (zones K, L, M and 
N). Since every two neighbor zones share one row or one column of photomultiplier 
tubes (PMT), one gamma hit may trigger two zones if its energy is collected by the 
shared PMTs. To find out the zone that has the "real" hit, a threshold may be set for the 
10 pre-normalized position signals of each zone to reject the "false" hit from its neighbor 
zone. 

FIG. 26 shows an example of a two-zone design. Each zone has one independent 
block detector coupled to 4 PMTs. The pre-normalized position of zone 1 is X1=A+D, 
and the pre-normalized position of zone 2 is X2^C+F. If a gamma hits zone 2, there is no 

15 contribution to XI . However, if a gamma hits zone 1 (even in the corner crystal located 
in the center of PMT-B or PMT-E), to have the crystal well decoded, it still has more than 
3-5% of total energy going to PMT-A and PMT-D. That means the pre-normalized 
position signal (A+D) will have more than 3-5% of total energy for a real hit. Therefore, 
in data acquisition with a threshold E0 for the total energy, if we set 3-5% of E0 for a pre- 

20 normalized position signal threshold, the false hit caused by its neighbor's zone will be 
rejected (since it does not meet the signal threshold amount). Using a threshold 
corresponding to the energy ! s threshold for pre-normalized position is an effective way to 
reject false events within a zone. 

The HYPER method combined with multiple zones may significantly increase the 
25 count-rate capability and dose efficiency of gamma cameras. For coincidence detection 
applications such as, but not limited to, positron coincidence imaging, a special additional 
technique is provided in this disclosure. Conventionally, integration time is fixed; 
therefore, the end of integration and signal-digitization (for energy and position) is 
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synchronized with the trigger signal (just delayed by the fixed integration time). The 
trigger signal is generated at the leading edge or the arrival time of the event. 
Coincidence detection can be performed between the ends of signal digitization or the 
trigger signals. But, HYPER uses a dynamic integrating method, which means the 
5 measured energy and position signals are no longer synchronized with the leading-edge 
of the present triggering signal because the energy and position signals are only generated 
or digitized at the arrival of the next event (arriving at a random time from present event). 

To detect a coincidence event using the HYPER method, a trigger signal may be 
delayed by a fixed time and a new synchronization may be setup between the delayed 

10 trigger and the energy/position signals before they can go to a coincidence procedure. 
The energy and position signals of all events may be buffered into a FIFO (first-in-first- 
out) queue. Then, a delayed trigger signal may be used to read them out individually to 
restore the synchronization between the events output of the FIFO. This process is shown 
generally in FIG. 27. The fixed delay time should be longer than the maximum 

1 5 integrating time. The trigger delay circuit should have good timing-resolution and a 
small dead time. 

The HYPER methodology of this disclosure, the multiple HYPER zones method, 
and the event-synchronization method for coincidence detection may all improve PET 
camera imaging performance and facilitate gamma cameras for radionuclide therapy 
20 dosimetry imaging, cardiac first-pass imaging, positron coincidence imaging, and the 

simultaneous acquisition of transmission and emission data using difference isotopes with 
less-contamination between transmission and emission data. Those having skill in the 
art, with the benefit of this disclosure, will recognize additional suitable uses for these 
techniques. 



25 BRIEF DESCRIPTION OF THE DRAWINGS 

The following drawings form part of the present specification and are included to 
further demonstrate certain aspects of the present invention. The invention may be better 
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understood by reference to one or more of these drawings in combination with the 
detailed description of specific embodiments presented herein. 

FIG. 1 is a graphical representation of scintillator light output of gamma rays over 

time. 

5 FIG. 2 is a graphical representation of scintillator light output of a plurality of 

radiation particles having pile-up. 

FIG. 3 A is a graphical representation of exponential decay of a radiation particle 
over time. 

FIG. 3B is a graphical representation of total integrated energy from a radiation 
10 particle. 

FIG. 4A is a graphical representation of a weighted instantaneous and integrated 
signal of a radiation particle according to the present invention. 

FIG. 4B is a graphical representation of a weighted sum of an instantaneous 
signal and an integrated signal of a radiation particle according to the present invention. 

1 5 FIG. 5A is a block diagram of an exemplary embodiment of the present 

invention. 

FIG. 5B is a block diagram of another exemplary embodiment of the present 
invention. 

FIG. 6 is a block diagram of an exemplary apparatus according to the present 
20 invention. 

FIG. 7 is an illustration of a PMT configuration that may be used in accordance 
with the present invention, 

FIG. 8 is a representation of boundary artifact generation. 
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FIG* 9 is a block diagram of an exemplary apparatus according to the present 
invention. 

FIG. 10 is an illustration of study results according to the present invention, 

FIG. 11 is a graphical representation of scintillator light output of a plurality of 
5 radiation particles having continuous pile-ups. 

FIGS, 12A-12B are graphical comparisons of energy spectra of count-rates of 
prior art methods (FIG. 12A) and methods according to the present invention (FIG. 12B). 

FIGS. 13A-D are graphical representations of energy spectrum of an apparatus 
according to the present invention at very high count-rates. 

10 FIG. 14 is a graphical representation of energy resolution (percentage error) as a 

function of count-rates of prior art methods and methods according to the present 
invention. 

FIG. 15 is a graphical representation of detected count-rates as a function of true 
count-rates of prior art methods and methods according to the present invention. 

15 FIG. 16 is a graphical representation of energy spectrum of a delay line pulse 

clipping method. 

FIGS. 17A-F are graphical representations of energy spectrum of a 99ra Tc source 
with and without remnant subtraction. 

FIG. 18A is a representation of a Monte Carlo result according to the present 
20 invention. 

FIG. 18B is a representation of a Monte Carlo result of a pulse clipping method. 

FIG. 18C is a representation of a Monte Carlo result of a conventional fixed 
integration method. 
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FIG, 19A is a flow chart of an exemplary method according to the present 
invention, 

FIG. 19B is a flow chart of another exemplary method according to the present 
invention. 

5 FIG. 20 illustrates an exemplary circuit according to one embodiment of the 

present disclosure. 

FIG. 21 illustrates a pileup condition according to one embodiment of the present 
disclosure. 

FIGS. 22-24 illustrate various calculation circuits according to embodiments of 
1 0 the present disclosure. 

FIG. 25 illustrates a four-zone detector according to an embodiment of the 
present disclosure. 

FIG. 26 illustrates zone identification techniques according to an embodiment of 
the present disclosure. 

15 FIG. 27 illustrates a coincidence calculation circuit according to an embodiment 

of the present disclosure. 

FIG. 28 illustrates a high resolution time lapse measurement circuit according to 
an embodiment of the present disclosure. 

DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 

20 The present invention in a broad aspect comprises a dynamic pile-up prevention 

technique for increasing the count-rate capability of scintillation detectors. In one 
embodiment, a Nal(Tl) scintillator may be used to illustrate the method. Although 
Nal(Tl) scintillators are discussed herein, the present invention is applicable to other 
detectors, such as BGO, GSO, LSO, plastic, and CsL The dynamic pile-up technique of 
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the present invention has a variable detector signal collection time (deadtime), whereas 
conventional systems typically have a fixed deadtime of 1 jus (4%) for Nal(Tl) 
scintillators. The reason that 1 |is is generally used as a fixed deadtime is that the 
scintillation light from a Nal(Tl) scintillator decays exponentially with a time constant of 
5 0,24 |o,s. Thus, approximately 98% of the light is collected in a 1 jis signal-collection 
time. Hence, energy measurement error (energy resolution) is minimized. 

From Poisson statistics, if a system is counting at a rate of R = 10 6 counts per 
second (cps) and if the signal integration (collection) time is 1 p,s, the fraction of 
non-pile-up counts are only e~ RT = 37%. Hence, approximately 2/3 of the counts are pile- 

10 up events. Though shortening the signal collection time below 1 |is would improve the 
fraction of nonpile-up counts, it would also cause two problems. First, shortening the 
time could increase energy resolution or error (good or small energy resolution is needed 
to reject scatter noise), as energy resolution is inversely proportional to the square root of 
light collection. Second, even if the collection time is shortened electronically, 

15 scintillation light continues to be emitted by the detector from previously detected gamma 
ray particles. These remnant signals would be added erroneously to the signal of the 
following gamma ray particle. 

Position sensitive detectors with Anger decoding algorithms are used in SPECT 
cameras and PET cameras (NaI(Tl), BGO and LSO systems) to reduce production costs. 

20 One drawback of Anger detectors is lower detection rates because all the photomultipliers 
(PMT's) involved in the localization of an incident event will be engaged in signal 
collection for a fixed time period, thereby inhibiting the detection system to process a 
second event incident within this time period. From energy and spatial resolution 
considerations, a fixed signal-collection time period of 2x s to 3x s is preferred, where t s is 

25 the scintillation decay-time constant. Due to the random distribution of the time lapse 
between two events, the average time lapse between two events should be 10 times the 
signal collection time (20x s - 30x s ) to reduce the probability of signal pile-up to 10%. For 
Nal(Tl) and BGO detectors, this limiting average time lapse would be 5-9 (as, which 
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corresponds to a maximum count rate of 110-200 thousand-counts/sec (Kcps) per crystal 
or crystal-block, which has been the maximum count-rates of most gamma-cameras and 
BGO block detectors for many years. However, for certain applications, higher operating 
count-rates are preferred. 

5 Recently, with the advent of coincidence positron imaging gamma-cameras, the 

need for increasing camera count rates has become more immediate. Currently, two 
approaches are used to increase gamma-camera count rates: 

(a) a dynamic integration method which integrates the first signal until a next 
event (pile-up) arrives. An estimated signal supplement is then added to the first pulse to 
10 correct for the signal collection deficit, and the same correction signal is subtracted from 
the next pulse to correct for the remnant light from the present pulse (Lewellen et aL 9 
1989). This method is relatively simple for a two event pile-up, but is more complicated 
for a multiple-event pile-up, where an additional circuit is needed for each higher 
multitude of pile-up. 

15 (b) a delay-line pulse-clipping (DLPC) method, which has been applied to 

Nal(Tl) PET cameras (Karp et aL, 1986). This technique has been adapted recently to 
coincidence imaging gamma-cameras. Generally, the DLPC technique reduces the 
scintillation pulse width to x (from 3x), which increases the maximum camera count rates 
by about 3 times. An advantage of this method is that it can be integrated into a 

20 multi-zone architecture, which splits a Nal(Tl) crystal into multiple pseudo-independent 
zones (as if having multiple pseudo-independent crystals) for increasing the maximum 
camera count rates further. It is known that the DLPC technique (Muehllenhner and 
Karp, 1986; Karp et aL, 1986; Tanaka et aL, 1979; Miyaoka et al, 1996b) may be used to 
alleviate the remnant signal problem, but such a technique does not remedy the degraded 

25 energy resolution problem. Furthermore, this technique only defers the onset of pile-up 
to a moderately high count rate. 

Another proposed approach disclosed in U.S. Patent No. 5,430,406, which is 
hereby incorporated by reference, uses a dynamic approach to obtain a weighted sum of 
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the instantaneous signal and the integrated signal to depict the energy for the first event in 
a two event pile-up. The second event in this two-event pile-up was only used for 
detecting a pile-up condition to process the first event, and this second event is not 
processed further for recovery attempts. Hence, this method may maintain a good spectra 
at moderately high rates when the most of the pile-ups are two-event pile-ups, but the 
count-loss will be high due to the loss of the second event. Second, it can be shown 
theoretically that the weighted sum is generally not equal to the energy of an impinging 
event (See Appendix), because it also includes the remnant light from all previous events. 
The weighted sum is only equal to the event energy when the first event in the two-event 
pile-up is not a pile-up on previous events, which is a condition in contradiction with very 
high count-rate situations, where most events are riding on the signal of one or more 
previous events. Hence, the errors in this energy measurement method will be large in 
very high count-rate situations (when many events are part of multiple-event pile-ups), as 
shown in the following table derived from Poisson's statistics for NaI(Tl) for the 
probabilities of multiple-event pile-ups within 1 us: 



TABLE 2 

PILE-UP FRACTION AS A FUNCTION OF COUNT-RATES 





no pile-up 


2-event 


3 -event 


4-event 


5 -event 


6-event 


7-event 






pile-up 


pile-up 


pile-up 


pile-up 


pile-up 


pile-up 


1 Mcps 


0.37 


0.37 


0.18 


0.06 


0.02 


0 


0 


2 Mcps 


0.14 


0.27 


0.27 


0.18 


0.09 


0.04 


0.01 


4 Mcps 


0.02 


0.07 


0.15 


0.2 


0.2 


0.6 


0.1 



The present invention provides a new approach for preventing signal pile-up, 
namely, a hybrid signal processing method based on several concepts. It is known that 
after a scintillation crystal detects a gamma ray, the light output decays exponentially. 
This signal is illustrated in FIG. 3 A, together with its one-standard-deviation (±a) error 
boundaries. Initially, about 40 electrons per 10 ns period are generated from a 140 KeV 
energy deposition (from 99m Tc). The total integrated signal over a period of 1 us is 
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approximately 1000 electrons from this 140 KeV energy deposition event, as shown in 
FIG. 3B. For a scattered gamma ray of 70 KeV energy deposition, the initial signal 
would be 20 electrons per 10 ns period, and the total integrated signal would be 500 
electrons. Although these numbers are specific for gamma rays detected by a Nal(Tl) 
scintillation detector, the same general principles, and the general curves shown in 
FIGS. 3A and 3B apply to other radiation particles and other radiation detectors. Thus, 
both the instantaneous signal and the integrated signal may be used to determine the 
detected gamma energy. However, the percentage-error is much less with the integrated 
signal collected over 1 us, as shown in FIG. 3B. For this reason, only the integrated 
signal technique is used traditionally. 

The present invention includes a hybrid signal processing technique that uses both 
the instantaneous signal and the integrated signal together to help derive the radiation 
energy, since both signals contain radiation energy information. 

In an exemplary embodiment, the instantaneous scintillation signal from the 
photomultiplier-tube (PMT) may be smoothed with a filter, for example a 10 ns 
time-averaging RC-filter. This smoothed instantaneous signal and the integrated signal 
are shown in FIGS. 3A and 3B. The smoothed instantaneous signal is amplified to the 
equivalent of 1000 electrons at time zero, and this amplified instantaneous signal is 
summed with the integrated signal (unamplified) as shown in FIG. 4A. Although the 
particular embodiment described herein amplifies this signal to a certain level of 1000 
electrons, it is to be understood that the smoothed instantaneous signal may be amplified 
to various levels, depending upon the type of radiation and detector. In certain 
embodiments, the smoothed instantaneous signal need not even be amplified. 

The resultant weighted-sum signal of the instantaneous signal and the integrated 
signal is a constant and is always a measure of the radiation energy if the event is not a 
pile-up on previous events, regardless of when the sum-signal is sampled. In an 
exemplary embodiment, the sampled signal amplitude may always be 1000 electrons, 
which is equivalent to a 140 KeV energy deposit. The weighted sum is constant in time, 
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because the amplified (weighted) instantaneous signal decays as T(E 0 e" t/ 7T) = E^ 1 (a 
property of Nal(Tl) and other scintillators), and the integrated signal increases as 
E o(l - the mathematical property of integrating an exponential function. Thus, the 
weighted-sum is independent of the time at which the signal is sampled by the 
components of the present invention. Combining the instantaneous signal and the 
integrated signal derives the weighted sum according to the following equation. 

Weighted Sum = E 0 e n + E 0 (l - e^) = E c (1) 

Given this equation for the weighted sum signal, it can be seen that at certain 
times, more of the signal will be derived from the instantaneous signal, while at other 
times, more of the signal will be derived from the exponential signal. If the weighted 
sum is sampled at earlier times, the variance of the measurement is higher, because more 
information weight is given by the instantaneous signal, as shown in FIGS. 4A-B. If the 
sum signal is sampled at 1 us, the sum converges to that of the integrated signal with a 
smaller variance, as shown in FIG. 4B. 

In the present invention, if there is no next gamma ray detected within the 
conventional integration time period, for example, 1 us, the sum signal acquired in the 
duration of the present event is sampled. If a next gamma ray arrives within 1 (as, the 
signal (sum) in the duration of the present gamma ray is sampled immediately before the 
scintillation signal of the next gamma ray reaches the signal measurement circuit. Hence, 
the measured signal sum relating to the present gamma ray excludes the erroneous pile-up 
signal from the next gamma ray. However, the signal sum may still contain signals from 
previous events. 

A flow chart of an algorithm according to the present invention is shown at FIG. 
19A. As shown in FIG. 19A, a detected signal is integrated in step 210, is amplified in 
step 220, and causes a trigger signal in step 230. The integrated value of the detected 
signal and the amplified value of the detected signal are summed in step 240. The 
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summed value is stored in step 250 and is sent also to the residual subtraction circuit for 
use in processing. 

In addition, the trigger signal from step 230 is sent to determine time lapse 
between the present detected signal and the previous detected signal in step 260, The 
5 sum of the previous detected signal is multiplied by an exponential of the time lapse in 
step 280. This value is then subtracted from the present detected signal sum value in step 
270. The output of step 270 is the detected energy signal E i? as shown in step 290. 

A block diagram of an exemplary embodiment of the present invention that 
implements this technique is shown in FIG. 5 A. As shown in FIG. 5 A, the exemplary 

10 pile-up prevention circuit 100 receives electrical signals from a photosensor (not shown 
in FIG. 5 A) connected to the input of the circuit. The pile-up prevention circuit (PPC) 
100 includes a timing trigger 110, which provides a signal to the sample and discharge 
control circuit 170 and residual subtraction circuit 190 upon receipt of a signal. In an 
exemplary embodiment the timing trigger 110 may be a simple threshold discriminator or 

15 a delay-line clipping discriminator, which minimizes retriggering by the remnant signal 
of the same particle. Alternately, the timing trigger 110 may be another appropriate 
triggering device, for example, a Schmitt trigger. PPC 100 also includes a smoothing 
circuit 120 to smooth the received signal. In an exemplary embodiment, the smoothing 
circuit 120 may be a filter, such as a low-pass RC filter or a slow amplifier. Further, in 

20 certain embodiments, a smoothing circuit may not be necessary. The smoothed signal 
passes to a delay circuit 130, which provides a desired delay to the signal before 
application to the calculation portion of the circuit. In an exemplary embodiment, the 
delay circuit 130 may be a delay cable or analog delay line integrated circuit. The delay 
circuit 130 may provide for a delay of between about 1 to 100 ns, and more preferably 

25 between about 5 to 50 ns. 

As shown in FIG. 5A, the calculation portion of PPC 100 includes amplifier 140, 
integrator 150, and adder 160. The amplifier 140 amplifies the smoothed signal, thereby 
providing a weighting factor to the signal. As discussed above, in an exemplary 
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embodiment the amplifier 140 amplifies the signal to 1000 electrons. In appropriate 
embodiments, an amplifier may not be necessary. The integrator 150 performs an 
integration of the signal over the exponential decay of the signal. The outputs of the 
amplifier 140 and integrator 150 are input into the adder circuit 160, which develops the 
5 sum of these signals. The resultant sum is forwarded to sample and discharge control 
circuit 170. 

In an exemplary embodiment, the sample and discharge control circuit 170 may 
be comprised of a register or the like, and a control circuit, such as a fast analog switch, 
FET transistor, or the like. Upon a signal by time trigger circuit 110, the sample and 
10 discharge control circuit 170 passes the resultant sum to analog to digital converter 

(ADC) 180. The ADC 180 converts the sum to a digital signal, which is then forwarded 
as the weighted-sum signal to the residual subtraction circuit 190, which subtracts any 
remnant signal from previous radiation events. 

As shown in FIG. 5 A, the residual subtraction circuit 190 includes a lapse timer 
15 192, which determines a AT (or time difference) corresponding to the time between 

receipt of two detected signals. The residual subtraction circuit 190 also includes a look- 
up table 194 that determines an exponential value based on the time difference. This 
exponential value is then multiplied by the previous sum signal S n (not shown). Finally, 
a subtractor 196 subtracts this remnant sum signal from the sum signal S { relating to the 
20 present detected signal to thereby output a digital signal E { corresponding to energy value 
of the present signal. 

After sampling, the sample and discharge control circuit 170 discharges the output 
of the integrator 150 to zero immediately by sending a control signal upon receipt of the 
signal from the timing trigger 110 to a capacitor and switch (not shown) associated with 
25 the integrator 150. 

In operation, the timing trigger 110 senses the arrival of the initial burst of 
scintillation electrons from any gamma ray. Immediately, time trigger 110 sends a signal 
to the sample-and-discharge control circuit 170 to sample the weighted sum of the 
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previous gamma ray and pass the sum signal to the ADC 180. Since the timing trigger 
110 has no delay, whereas the sum signal processing branch is delayed (in an exemplary 
embodiment, the delay may be between about 1 and 100 ns, and more preferably be 
between about 5 and 50 ns) with a delay circuit 130, the weighted-sum signal is sampled 
5 before the arrival of the energy signal from the next gamma ray. The delay is used to 
ascertain that the weighted sum of the present event is sufficiently discharged before the 
next event enters the integrator 150. 

Hence, pile-up of the next signal onto the present is prevented, and the weighted 
sum will not contain the energy of the next gamma ray. The delay also provides enough 
10 time delay for the sample-and-discharge control circuit 170 to discharge (to zero) the 
previous integrated signal from the integrator 150 before the arrival of the next signal, 
which prepares the integrator 150 for the next signal. The discharge avoids pile-up of the 
old integrated signal onto the new integrated signal. 

Since it takes a finite time for the electrical-charge output of the integrator 150 to 
15 be fully discharged to zero, especially if the previous event has high energy, this 

discharging time is the ultimate practical limit (processing dead-time) restricting the 
maximum count-rate performance of the present invention. In an exemplary 
embodiment, as shown in FIG. 5B, a dual-integrator design ("ping-pong") may be used as 
the integrator, with each integrator 150A and 150B taking turns in processing each 
20 consecutive event as follows: when a present event is detected, the circuitry will be 
switched via switch 145 to integrator 150A to integrate the signal of the present event; 
when the next event is detected, integrator 150 A will stop integrating the present event 
and the integrator output will be sampled. At the same time, the circuits will switch the 
second signal to integrator 150B for integration without waiting for integrator 150A to be 
25 discharged. Hence, integrator 150 A is discharged in this idle duty cycle when the 

integration of the next event is performed by integrator 150B. This will provide more 
than enough time for the integrator 150A to discharge fully. 



-23- 

25075668.1 



When a third event is detected, the signal will be switched back to integrator 
150A, which is fully discharged, for integration, and the integrator 150B is sampled and 
discharged at its leisure. Hence, with a dual-integrator "ping-pong" design, the 
discharging dead time is eliminated (to increase processing count-rate) and the error in 
energy measurement due to incomplete discharge is also eliminated. 

As shown in FIG. 5B, there may be two independent channels, each of which may 
comprise an integrator 150A or 150B, an adder 160A or 160B, a sample-and-discharge 
circuit 170A or 170B, or an ADC 180A and 180B, together with an input switch 145 and 
an output switch 185 for selecting the processing-channel. Each channel takes turn to 
process the next incoming event. In other embodiments, there may be more than 2 
channels present. This also provides an idle duty cycle for the integrators in each channel 
to be sampled and fully discharged. This dual-channel design also allows the use of a 
slower ADC (50% slower) to digitize the signals thus lowering cost. 

The digital remnant subtraction circuit 190 may include a digital look-up-table 
(LUT) for determining the remnant scintillation energy (from all preceding gamma rays) 
that the scintillator is still emitting after the arrival of the latest gamma ray. Specifically, 
the remnant subtraction circuit 190 may include a register (not shown) to store a previous 
weighted sum, a LUT 192 for generating an exponential term, and a digital multiplication 
IC (not shown). The input to the LUT is the trigger-time of the present and the preceding 
gamma ray (actually, their time difference At). The digital energy output of the remnant 
subtraction circuit 190 is thus the pile-up free energy (Ej) of the detected gamma, 

E-Si-S,^ (2) 

where S, is the weighted-sum relating to the present gamma ray, Sj_, is the weighted sum 
relating to the preceding gamma ray, and At is the time interval between arrival of the 
present and the preceding gamma ray. S; and S M include the remnant signals from all 
preceding gamma rays. This is an exact equation, regardless of how many multiple pile- 
ups exist in each sum signal (see APPENDIX for the mathematical justification). The 
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correction term and the exactness of this algorithm, being independent of the number of 
multiple-pile-up remnants still being emitted by the scintillator, are significant 
consequences of the present invention. Thus, the scintillation detector used in connecting 
with the present invention permit count-rates that produce 60% multiple pile-ups, since 
5 all the remnants from multiple pile-ups can be easily corrected for by equation (2). The 
present invention processes all incoming events the same way, regardless of whether it is 
a pile-up onto a previous event or whether there is a pile-up from the next event. 

All the components in a circuit according to the present invention are basic 
required components of conventional detector electronic systems, with the exception of 

10 the timing trigger and the residual subtraction circuit. Therefore, the real-time pile-up 
prevention/correction provided by the present invention is an inexpensive solution to the 
problems inherent in the prior art. If there is no pile-up, the system converges to a 
conventional integrating circuit, as the sample and discharge control circuit 170 may be 
set to pass the resultant sum signal at a time approximately equal to 4t (i.e., 

15 approximately 1 jus for the Nal(Tl) detector discussed herein). In case of a pile-up, it 
correctly measures the energy of all the pile-up gamma rays. 

As long as the average time lapse (AT) of arrivals between two gamma rays is 
longer than the average detector integration time (Tc = 1 jas), the average fraction of 
scintillation light collected will not diminish, which means the detected energy resolution 

20 will not be compromised. Hence, a Nal(Tl) scintillation detector employing the methods 
and apparatus of the present invention has the potential to count up to 1 x 10 6 cps without 
degrading the average energy resolution. A traditional detector circuit with a deadtime of 
1 ills must lower the acceptable count-rate by 10 times, to 100,000 cps, to prevent the 
random arrival of the next gamma ray within 1 jus of the present gamma ray (Nicholson, 

25 1974). Poisson statistics shows that the non-pile-up fraction is 90% for 100,000 cps, 
which is why conventional Nal(Tl) detection systems cannot count over 100,000 cps. 
The present invention can count even faster than 1 x 10 6 cps with some compromises in 
energy resolution. 
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The PPC of the present invention may be adapted for use with many different 
detectors and imaging systems. For example, the present invention may be used with 
thyroid probe, position sensitive detectors, or gamma ray detectors which detect neutrons, 
charged particles and gamma rays. 

The conventional integration method (and its equivalent pulse-shaping method) 
has been the gold standard in the nuclear physics and nuclear medicine field for the past 
40-50 years. The present invention is a significant breakthrough. This breakthrough can 
be applied to many nuclear detection areas, such as nuclear physics and engineering, high 
energy physics, nuclear medicine, industrial gauging, and oil field down-hole logging, 
among others. Certain applications will be discussed presently. 

Nuclear Medicine Applications 

The present invention has special attributes that make it particularly applicable to 
medical imaging, and especially gamma cameras. When a gamma ray hits a detector in a 
gamma camera, both energy and position information are needed for creating the image. 
The energy information discriminates (or qualifies) the useful imaging gamma ray from 
background gamma rays which have been scattered around in the body being imaged 
(scattering lowers the gamma energy). When a detected gamma is qualified by its 
detected energy, the position information locating the point of detection of the gamma ray 
has to be determined to create a gamma ray map (image). 

The pile-up prevention circuit (PPC) of the present invention can determine the 
gamma ray energy, despite pile-ups. However, the present invention would be less useful 
for gamma cameras if it could not also determine the gamma ray position in pile-up 
situations. A particular embodiment of an apparatus according to the present method may 
be incorporated into a regular gamma camera to determine position information, as well 
as energy information in a pile-up situation. 

This embodiment may be known as a position-energy pile-up prevention (PEPP) 
algorithm and circuit. This solution is compatible with existing gamma camera 
electronics. Compatibility is important because it allows the pile-up-prevention circuit of 
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the present invention to be applied to present gamma cameras without major modification 
to their existing front-end electronics. 

A regular gamma camera and its electronics generate five signals of interest for 
the present invention, namely X + , X_, Y+, Y_ 5 and Z. Z is the energy of the gamma ray. 

The X-position and Y-position of the gamma ray can be calculated from these five 
signals, 

X = (X + -X )/Z (3) 
and 

Y = (Y + -Y)/Z (4) 

The prenormalized positions are defined as X' = (X + - X_), and Y* = (Y + - Y_). Hence, the 
X-position and Y-position can be calculated by a subsequent normalization with Z: 

X = X'/Z (5) 

and 

Y = YVZ (6 ) 

hi addition to preventing energy pile-ups in scintillator-PMT's, the same PPC 
circuit may be applied directly to the prenormalized position signals, X* and Y', to prevent 
pile-ups by generating pile-up-free X' and Y'. The mathematical proof of this property 
can be found in the APPENDIX. With the realization of this important property, the PPC 
circuit of the present invention may be connected to the prenormalized positioning 
signals X' and Y' in a regular camera, as well as the energy signal Z' (a fast instantaneous 
signal sum of all the PMT's), as shown in FIG. 6, to generate three pile-up-free signals 
X', Y' , and Z' for calculating the energy and position of the detected gamma ray despite 
multiple pile-ups. 
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The system of FIG. 6 includes a conventional gamma camera 200, which develops 
prenormalized position signals X' and Y', and energy signal Z'. Each of these signals is 
input into a separate pile-up prevention circuit (PPC) 210. This PPC 210 may be the 
same circuit as discussed above, in connection with FIGS. 5 A or 5B, or it may be a 
5 slightly modified version (as shown in FIG. 6, certain components of the circuit of FIG. 
5A have been deleted for illustrative purposes). First the PPC 210 may not have a time 
trigger, as trigger signals for this embodiment may come directly from a fast trigger 220, 
which is input to all of the PPC's 210. This fast trigger 220 is a new addition to a 
conventional gamma camera. Further, the individual PPC's 210 may not have a remnant 

1 0 subtraction circuit, as all signals output from each PPC 210 may be input into a digital 
signal processor (DSP) 240. In an exemplary embodiment, the DSP 240 may include the 
remnant subtraction algorithm for the X', Y\ and Z signals. The DSP 240 provides for 
faster processing. However, in certain embodiments, it may be possible to have the 
lookup tables within the individual PPC 210 units. The DSP 240 may also include other 

1 5 corrections (energy and distortion corrections) for the X', Y' and Z signals, and X, Y 
renormalization process. 

The system of FIG. 6 also includes a coincidence timer 230, which receives a 
trigger signal from the fast trigger 220. Upon receipt of the signal, the coincidence timer 
230 sends a time clock signal to FIFO's 235, which is transmitted to the DSP 240. This 
20 time clock signal is for pile-up remnant subtraction and for coincident timing 
measurement in positron coincidence imaging (PET). 

As discussed above in connection with FIG. 5 A, the prenormalized position 
signals (X', Y') obtained from gamma camera 200 are delayed, amplified, and integrated 
within the PPC 210. The amplified instantaneous position signals xX'(t) and xY'(t) are 
25 summed to their integrated position-signal, 1 X' (t)dt, \ Y' (t)dt, respectively, also within 
the PPC 210. These weighted sums of prenormalized positions are immediately sampled 
by a fast analog-to-digital converter (ADC) within PPC 210, just before the onset of the 
next pile-up event. The arrival of the next gamma ray is sensed by the fast trigger 220, 
which monitors the fast instantaneous total energy output, Z'. After the ADC sampling, 

-28- 

25075668.1 



the integrator circuit of PPC 210 is immediately discharged so that it can start integrating 
the position signals coming from the next pile-up gamma ray. Alternatively, there may 
be two matching integrators taking turns in processing each consecutive signal to allow 
for an idle duty cycle for discharging the integration to zero. Further, as discussed above, 
5 there may be two matching channels of amplifiers, integrators, samplers, and ADC's to 
handle each successive event to allow the integrators in each channel to discharge to zero. 

The fast trigger 220 also marks the arrival time of all gamma rays. The timing 
mark is generated by a fast clock, which in an exemplary embodiment may be a 250 MHz 
clock. These arrival times are sent to the DSP 240 for remnant signal subtraction and for 
10 positron coincidence detection in PET imaging. The remnants of position signals, 
generated by the residual light output of all preceding pile-up gamma rays, will be 
subtracted from the present position weighted-sum in the DSP 240 in FIG. 6 using a LUT 
that stores the exponential term in the following remnant subtraction operation to 
determine the remnant position signal to be subtracted from they-th event: 

15 RP^SP^e^ (7) 

and 

Xj^SPj-RP, (8) 

where S? jA is the last position-weighted sum and At is the time-lapse between the j and 
the j-1 position signals. This remnant subtraction strategy automatically subtracts for all 

20 higher levels of multiple pile-ups. After the remnant subtraction, the DSP 240 calculates 
the normalized positioning estimation, X = X' / Z and Y = Y7 Z , where Z is the total 
energy after remnant correction of the present event. Alternate methods of remnant 
subtraction may include performing the process within the individual PPC 210 units. The 
DSP 240 also performs the regular correction processes for linear distortion, field 

25 nonuniformity, and regional signal-pulse-height variation, in real time. The host 
processor 250 is for image display and operator interface. 
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An apparatus according to the present invention thus decodes the position and 
energy of each detected event. The fast trigger 220 detects the arrival of an event, and 
triggers a time-mark output in the coincidence-timer 230 for the event's arrival time. This 
event-arrival time mark will be used to calculate the prenormalized positions and energy. 
5 The calculation may be performed by the DSP 240, where exponential functions may be 
stored as a look-up-table. The DSP 240 also performs the energy normalization on the 
prenormalized positions to provide the true position signals. In this processing scheme, 
all the events (photopeak or true events, and scatter noise events) are included in the 
processing to calculate all the remnant signals from previous events, and an event's real 

10 energy is not known until after processing. Hence, the energy acceptance of an event has 
to be made after the decoding operation, in/after the DSP 240, If the decoded energy is 
higher than that of the scatter radiation noise signal, the detected radiation is accepted as a 
true (photopeak) event. The output of the DSP is thus the pile-up-free, 
energy-normalized positions (x, y) and gamma-energy (z). Since all data acquisition 

15 computers in gamma cameras accept the standard signals of (x, y, z), the proposed 
processing algorithm and electronic architecture is compatible with existing data 
acquisition computers. This compatibility is useful to adapt the scheme to existing 
cameras. 

Implementation of the present invention in a conventional camera is very feasible, 
20 as shown in FIG. 6. A DSP circuit is standard in most cameras for processing position 
signals and for distortion correction. Thus, the only additional circuits needed for 
incorporating the methods of the present invention into a regular camera are the three 
channels of PPC circuits according to the present invention, a fast trigger circuit, and a 
clock, which increase production costs by only a few hundred dollars to achieve a 10-20 
25 times increase in count-rate capability. 

Such an increase in count-rate capability would permit or improve the following 
imaging procedures: 
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(a) Positron imaging, because it requires the removal of a lead collimator 
from a gamma camera, which exposes the Nal(Tl) detector to a 10-fold or more 
gamma-ray flux. This forces an injected positron tracer dose to be reduced by 5-1 Ox to 
prevent significant pile-up in a regular camera. With the present invention, it would not 

5 be necessary to reduce the injected dose. The positron image quality would be 

significantly improved as a result of this 5-10 times increase in positron counts in the 
image. Imaging time can also be decreased as a result of the increased count-rate 
capability; 

(b) Imaging radionuclide therapy patients to deduce actual radiation dose 
10 delivered to tumors and organs so that radionuclide therapy can be improved; 

(c) Dynamic first-pass cardiac imaging to study shunts, valves, the right heart 
and lung; 

(d) The use of very short half-life tracers to reduce radiation dose and increase 
patient throughput; 

1 5 (e) The use of larger Nal(Tl) detectors to reduce scanning time for 

whole-body imaging; 

(f) Acquiring very high count-rate transmission data while the emission data 
is being collected. In such a scheme, a lower energy gamma source is used for the 
transmission data; the very high-rate transmission collection is used to minimize the 
20 effect of down scatter 1 contamination from the emission gamma. Other improvements in 
medical imaging would also be produced by application of the present invention. 

The present invention also provides for further improvement of the detection rates 
of a gamma camera. In this embodiment, the camera signals may be split geometrically 
into 4 or more independent signal processing zones as if four independent cameras 
25 existed adjacent to each other. More specific details of such a system are disclosed in 
U.S. Patent No. 5,319,204 and, U.S. Patent No. 5,453,623, the disclosures of which are 
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hereby incorporated by reference. The incorporation of the PEPP (position and energy 
pile-up prevention) algorithm of the present invention into such a multi-zone design 
would further increase the count-rate capability of the camera by two or three times. 

Multi-zone designs have been proposed and implemented in NaI(Tl)-based PET 
cameras to increase maximum count-rate capability (Muehllehner et al, 1995; 
Muehllenhner and Karp, 1986; Karp et al, 1986; Freifelder et al y 1994). These designs 
depend on the fact that most of the scintillation light is distributed only to the neighboring 
7-9 PMTs. 

For a camera with 36 PMT's (6x6 configuration), the camera may be divided 
into four identical square zones, identified as K, L, M, N, as shown in FIG. 7. Each zone 
therefore has nine PMT's (3 x 3). Each zone may be treated as an independent camera 
with its own PEPP circuit according to the present invention to correct for pile-ups and 
image distortions. The nine PMT's in each zone may be grouped into one 
signal-triggering line. With four zones, there would be four signal-triggering lines, S , 

K. 

S l' S m' S n If the first S amma ra y is detected in zone K, the PEPP circuit of zone K 
would be turned on to measure its position and energy. If a second gamma is detected in 
zone K within the pile-up time, the PEPP circuit would correctly measure the position 
and energy of both gamma rays. 

If the second gamma ray strikes zone L while zone K is processing a prior count, 
two scenarios may occur. These two scenarios are addressed in detail: 

(1) In the first scenario, the gamma rays detected in zone K and zone L are 
both far from the zone boundaries between K and L, such that there is little light spilled 
from one zone to the next. Hence, there is little or no signal interaction between the two 
events to cause positioning errors. Both gamma rays may be processed independently by 
the PEPP circuit of each zone to obtain the correct energy and position. 

(2) In the second scenario, the first gamma ray detected in K is so near the 
K-L boundary that the PEPP circuits in both zone K and zone L are triggered, as shown in 
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FIG. 8. In this case, zone K may decode a valid position, but zone L would decode an 
artifact-count at the K-L boundary near the event in K, as shown in FIG. 8. When the 
second gamma hits zone L later, at location A in FIG. 8, the remnant signal of the 
fictitious count created in the L-PEPP circuit would be subtracted from the signal of the 
5 second gamma ray (A), so that there is no interference in determining the position of the 
second gamma ray. Thus, in this scenario, both the first (K) and second (L) gamma rays 
would be correctly measured, but an additional artifact count will be created at the K-L 
boundary. 

With these inter-zone light-spill considerations, the processing circuitry for the 
10 multi-zone pile-up prevention (MZPP) camera is shown in FIG. 9. As shown in FIG. 9, 
each zone may have its own PEPP circuit 310. In addition, each zone may have its own 
DSP 315. These four PEPP circuits may be monitored by an inter-zone coincident 
detection circuit 330 and a multi-trigger processor 340. The interzone coincident 
detection circuit receives input from the fast timing triggers of each of the zone PEPP 
15 circuits 310. The multi-trigger processor 340 receives zone position signals from each of 
the zone PEPP circuits 310. In an exemplary embodiment, the processor 340 may be a 
DSP. The output of processor 340 is input into a main DSP processor 350, which 
provides information to image buffer 360, interface 370 and Host CPU 380. The main 
DSP processor 350 serves two functions: (1) it merges the four zone-images into one 
20 image, and (2) it performs a final distortion correction for the combined image. 

To eliminate the problem of additional artifact counts at a border, the interzone 
coincidence-triggering circuit 330 (between any two or more zones) may be used to 
detect the simultaneous triggering of two or more zones caused by an event detected near 
the boundary. When two or more neighboring zones are triggered by a single event near 
25 the boundary, the PEPP circuit 310 in each zone would respond (as originally designed) 
as if an independent event is detected in each zone by generating their own position 
signals, 

X =(X -X )/Z (8) 
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and 



X =(X -X )/Z (9) 

L L+ L- L V 7 

Furthermore, the simultaneous triggering of both zones would also activate the 
inter-zone coincidence circuit 330 which would then feed these two or more independent 
5 position signals into multi-zone-trigger processor 340, which performs a centroid 
averaging for these two or more position signals. Mathematically, this averaging is 
equivalent to combining the two or more zones and using the regular Anger positioning 
method over this larger domain, as shown in the following equation: 

X = (X K Z K+ X L Zj/(Z K+ Zj = {(X K+ -X K J + (X L+ -X L J}/(Z K+ Zj (10) 

10 This averaged position is stored in the main processor 350 that stores the 

composite camera image (four zones combined), whereas the individual zone position 
signals (X r and X ) are discarded. This method thereby eliminates the boundary-artifact 
event and also provides a better estimation of the position, since all of the scintillation 
light emitted would be used for computing the position (including light spilled into the 

15 adjacent zone). 

Inter-zone light spill also exists in the two high count-rate cameras using a 
multi-zone design because, as long as a single Nal(Tl) crystal is used, light will be 
distributed from one zone to the next. Technical solutions are readily achievable as 
indicated by AD AC MCD and UGM SPECT and PET cameras (Muehllehner et aL, 1995; 
20 Glass et aL, 1996; Freifelder et aL, 1994). The present invention, using a monitor to 
detect coincidence triggering of two zones is simple, straightforward to implement, and 
usable with the pile-up prevention technology of the present invention. 

An increase in detection rates by 2-3 times is expected over the single-zone 
implementation, if the distribution of the scintillation light is mainly limited to a circle 
25 4.5" in diameter centered on the point of detection. When this light-spread condition is 
met, only the gamma rays hitting the shaded area shown in FIG. 10 can trigger or affect 
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both zones. Hence, the probability that a gamma ray will hit the independent area is 
21/36, while the probability that it will hit the light-spill-over area is 15/36. Thus, the 
count-rate enhancement obtained by splitting the camera signals into 4 zones is (21/36) 
x4 or 2.3 times. 

5 Shown in FIG. 19B is a flow chart for the method of the present invention relating 

to this multi-zone pile-up prevention embodiment. As shown in FIG. 19B, an 
instantaneous detected signal is provided to four independent PEPP's in step 410. Then it 
is determined whether this signal creates a multi-trigger event in step 420. If it is not a 
multi-trigger event, control passes to step 430, in which it is determined whether it is a 
10 single-trigger event. If it is a single-trigger event, control passes to step 440 and single 
active zone position processing is performed, as discussed above. 

If it is determined that it is a multi-trigger event, in step 450 active zones are 
combined into a larger domain for processing. Then, the large domain position 
processing is performed in step 460, as discussed above. The result of either path of this 
15 method results in a determination of the true and accurate position of the instantaneous 
signal in step 470. 

Combining the pile-up prevention method and the multi-zone method, the 
count-rate of a camera according to the present invention may reach 
4,000,000-5,000,000 cps. With optimization and further development of the PEPP circuit 
20 as discussed earlier, the count-rate may be extended further. Monte Carlo simulation 
indicates that 2 x 10 6 cps is achievable with the pile-up prevention method alone, and a 
combined count-rate of 2.3 x 2 x 10 6 = 4.6 x 10 6 cps may be achievable, a significant 
improvement over the 10 5 cps achievable in a regular camera. 

The present application of scintillation detectors can be divided into four main 

25 areas: 

(i) industrial gauging, including thickness, level and density gauging; 
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(ii) surveying, including oil field logging, nuclear reactor monitoring, 
nuclear fuel cycle monitoring and airport luggage inspection for explosives; 

(iii) research applications in nuclear physics, high-energy physics, 
medicine, and industry; and 

5 (iv) medical uses (gamma cameras, PET cameras, bone scanners, 

thyroid probes, general purpose probes, monitors, and dosimeters). 

A new and potentially important application of the present invention is in airport 
surveillance to look for plastic explosives in luggage using, for example, neutron 
activation techniques. The ability to detect radiation 10-20 times faster would allow (a) 
10 detectors to be used in a much higher radiation areas to extend the usefulness of the 
detector system, (b) the ability to count faster would increase the speed of data 
acquisition to shorten the data collection time, e.g., in an airport bomb surveillance 
system, a 10-20x stronger neutron source can be used to scan luggage faster, so that 
luggage scanning times can be reduced by 10-20 times to improve throughput. 

15 In industrial gauging application, the scintillation detector is used in conjunction 

with radiation sources to measure the quantities or density of materials being processed. 
For example, while molten cement or molten metal is poured into a container, a nuclear 
level gauge may be used to stop the pouring process when a specific level is reached. 
Another example is "thickness gauges" using gamma-ray backscatter detectors in 

20 thin-film processing in the plastic and paper manufacturing industry (especially for 

high-value films such as video-tapes and electronic capacitor films). Quite often, the film 
to be scanned is manufactured in large-area sheets, and the detectors together with the 
radiation source must scan the entire sheet. The detector's ability to detect gamma rays at 
very high count-rate (10-20x) means that a very intense gamma-ray source (10-20x) can 

25 be used, which in turns translates to a 10-20x faster data acquisition. Thus, the time spent 
on the inspection-scanning processes or manufacture-control processes can be reduced by 
90%, thereby increasing production rate. 
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Scintillation detectors are also widely used in radiation monitoring in nuclear 
reactors, and in oil field borehole logging to survey rock/hydrocarbon structures along the 
borehole in oil exploration. A faster counting system would increase the information 
collected or decrease the data collection time, which would lower the cost of data 
acquisition (oil rig time is very expensive). 

Scintillation detectors are also widely used for nuclear applications, for example: 
(a) reactor monitoring of liquid and gaseous streams to look for isotopes, (b) fuel rod 
cladding failure, (c) isotope scanning for irradiation fuel to determine power distribution 
and migration of fission products in the reactor core, (d) reactor fuel fabrication and 
quality control, (e) spent fuel reprocessing, and (f) management of nuclear waste. Since 
these nuclear-reactor related detectors are used in very high radiation flux areas, this high 
count-rate invention is potentially very useful in the nuclear energy business. 

In research laboratories, a faster counting detector system is always welcome, as it 
allows new experiments to be performed. Nuclear-electronics instrumentation companies 
that market detector electronic modules for research laboratories may wish to add this 
electronic invention to their catalog. 

In medical imaging applications, present gamma cameras are adequate for 
imaging diagnostic quantities of single photon tracers, but their limited count-rate 
capabilities introduce image artifacts and count-loss, and degrade image quality when 
gamma ray flux is high. These situations may occur in; 

(a) positron imaging, which exposes the Nal(Tl) detector to a 
10-fold or more gamma-ray flux, with the removal of the lead collimator 
from a gamma camera for positron imaging; 

(b) imaging radionuclide therapy patients to deduce the actual 
radiation dose delivered to tumors and organs for improving radionuclide 
therapy treatment; 
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(c) dynamic first-pass cardiac imaging to study shunts, valves, 
the right heart and lung; 

(d) the use of very short half-life tracers to reduce radiation 
dose and increase patient throughput; 

5 (e) the use of larger Nal(Tl) detectors to reduce scanning time 

for whole-body imaging; and 

(f) acquiring very high count-rate transmission data while the 
emission data is being collected (the very high-rate transmission collection 
minimizes the effect of f down scatter' contamination from the emission 
gamma). 

With this invention, the count-rate capability of gamma cameras can be enhanced 
by 20-40 times, which would allow (or improve upon) the above useful imaging 
protocols. The new capability may reduce reliance on the use of PET costing $2,600,000 
to purchase and $250 3 000/yr to operate. Another application of the present invention is 
for the low-cost, high-resolution PET detection system of PMT-quadrant sharing design, 
as shown in U.S. Patent Nos. 5,319,204 and 5,453,623, to compensate for its lower count- 
rate. 

Thus, the present invention will open up new and exciting clinical applications for 
gamma cameras, which is especially important in today's healthcare environment when 
clinics must perform more with less money. 

HYPER Techniques and Applications 

As the role of gamma cameras expands in positron coincidence imaging, 
radionuclide therapy dosimetry imaging, and cardiac first-pass imaging, there is a need to 
significantly increase the count-rate capability of gamma cameras because of the large 
25 photon flux in these procedures. Conventional detector electronics and the traditional 
Anger-positioning algorithm, however, hinder higher count-rate imaging because of the 
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pileup of signals of the detected events in the detector. Multiple events may be merged 
into one artifact event in the detection system; hence, true events may be lost while 
artifact events are generated. The pileup problem also exists in dedicated PET cameras, 
albeit to a lesser extent. This disclosure allows one to recover a pileup-free pulse signal 
5 even in multiple-pileup conditions. 

This disclosure provides several ways to implement the formulae developed in the 
HYPER method, thereby either improving the energy and position resolution or 
increasing the processing speed for pile-up prevention techniques. This disclosure also 
introduces a zone-identification methods for the application of multiple HYPER circuits 
10 to further increase the count-rate in gamma and PET cameras. It also explains how to use 
the HYPER method in the application of coincidence imaging. 

The HYPER method involves a set of formulae that solve the signal pile-up 
problem in scintillation detectors for both energy and position signals despite the pileup 
of multiple signals. The spirit of the methodology is as follows: (a) dynamic integration 
15 - the signal integration triggered by a present-event stops immediately before the next 
event is detected; (b) calculating a weighted value for estimating a total energy in the 
scintillation, which includes the energy of the current event and a residual energy from all 
the previous events; and (c) remnant correction is used to calculate a pile-up free energy 
from two consecutive weighted values. 

20 The above three steps may be used to create a pile-up free signal. This disclosure 

explains several ways to implement and to improve the measured resolution of the 
weighted value for estimating the total energy by: (1) applying an analog filter to the 
weighted-sum of an instantaneous value and an integrated value; the weighted-sum signal 
is multiple sampled, and another digital filter may be applied to those samples to get a 

25 better weighted value; (2) dynamic, digital weighting the integrated value; (3) dynamic 
weighting the sum of the multiple-sampling (repeated sampling) of the values of the 
instantaneous signal. 
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In certain embodiments of this disclosure, a dynamic integrating method is used, 
which means that the measured energy and position signals of the present-event are no 
longer synchronized with the trigger. This is the case because the energy and position 
signals of the present-event are measured at the arrival-time of the next-event (at a 
5 random time), and the trigger is issued at the arrival time of the present-event. To detect 
a coincidence event using the HYPER method, the original trigger signal may be delayed, 
and a new synchronization may be setup between the trigger and the energy or position 
signals before they go to a coincidence procedure. The improved HYPER method, 
multiple HYPER zones method, and this synchronization method may all facilitate the 
10 use of PET for applications such as, but not limited to, positron imaging, gamma cameras 
for radionuclide therapy, dosimetry imaging, cardiac first-pass imaging, positron 
coincidence imaging, and the simultaneous acquisition of transmission and emission data 
using difference isotopes with less-contamination between transmission and emission 
data. 

15 Theory 

Every exponentially decreasing pulse signal, for example from a scintillation 
detector, includes a present pulse signal and a remnant signal from all the previous pulse 
signals. With reference to FIG. 2, the total light signal after event 4 is the sum of the light 
signal of event 4 and a remnant light signal. 

t r r 

where S n is the total energy inside the detector after event-n without post-pileups, E n is the 
energy of the event n, R^ is the remnant energy from all the previous events, and x is the 
decay time. So, we have S n = E n + R n . 
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By the time of t*.,, the pulse signal becomes -?-e~ ( '~'" )/T = s. e <^-0^ e <*-^Vr _ 

T T 

This pulse signal will become the remnant signal to event (n+1). That means 

?^L e -(t-i^ir = ^ Le -(w-0/^-(^„ +1 )/ r> Thus, we have R n+l = S n e {t »^' Vr . 

T T 

Therefore, we have S n+l = E n+1 + R n+l = E n+l + S n e~ {tn+X ~ t,l)lz . So, the energy of the event 
5 (n+1) can be calculated from the following equation: 

*. + i = 5. + i-V (, ' +, "' )/r (2-1) 



This method is good for every exponentially decreasing pulse signals with multiple 
pileups. It also applies to the pre-normalized position signals in Anger positioning gamma 
cameras, since the pre-normalized position signals are exponentially decreasing too. 

PX n+l =SX n+l -SX tt e^-^ (2.2) 
PY^^SY^-SY^'-*^" (2.3) 

10 

where PX^, PY n+1 are the pileup-free pre-normalized position signals for event (n+1) and 
SX^, SY a+1 are the total estimated pre-normalized position signals. 

Circuit design 

The present disclosure, in a broad aspect, comprises a dynamic pileup prevention 
1 5 technique for increasing the count-rate capability of Anger type cameras with scintillation 
detectors such as, but not limited to, Nal(TI), BGO, GSO, LSO, plastic, and Csl. The 
dynamic pileup prevention technique has a variable detector signal collection time 
(deadtime), whereas conventional systems typically have a fixed deadtime (4x, l|is for Nal 
scintillator). The integration of the first event stops immidiately if the next event is detected 
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with a certain period of time (4x). FIG. 20 shows an exemplary diagram of an appropriate 
circuit. 

Computation circuit 

The computation circuit in FIG. 20 calculates a total value of a first and second 
5 incoming pre-normalization signal and a third incoming energy signal. 




Weighted-sum signal 

Equation (2.4) indicates that the weighted-sum signal of the integrated signal and 
the instantaneous signal amplified by x is a constant. The weighted-sum signal is always a 

10 measure of the total radiation energy before the next event is detected, regardless of when 
the sum-signal is sampled. Since the instantanous signal is more noisy than the integrated 
signal, a low-pass smoothing circuit may be applied to the weighted-sum signal to reduce 
the noise. To further increase the signal to noise ratio, the weighted-sum signal may be 
multiple/continously sampled, and a digital filter may be applied to the samples to get a 

1 5 better estimated value. 

In one embodiment, the filter may be a median value filter or a weighting output of 
these samples. For example, S H = W X S\ + W 2 S\ + ... + W m S™ , where S l n indicates the i-th 
sample of the weighted-sum signal S n may be used. The sum of the weighting factor is (Wj 
+W 2 +. . .+W m =1). Since a later sample of the weighted-sum has better stastics than an 
20 earlier sample, the weighting factor for an earlier sample should be smaller than a later 

sample. For example, with W t = W i+l e~ T/T (where T is the time lapse between samplings), 
the weighting factor follows the same exponentially decay as the pulse signal. 
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FIG. 24 shows a diagram of an appropriate decay filter. The number of samples (N) 
is not fixed. Thus, to obtain the estimated total energy, the accumulated output signal may 
be normalized by: 

i -N*AT/t 



5 In one embodiment that leads to a a simple electronics design, one may buffer the 

several latest samples of the weighted-sum signals and create an average of the buffered 
signals for the estimated total energy. 

Weighted integrated-value 

Equation (2.5) indicates that the total estimated energy may be calculated by 
10 weighting the integrated value (integrating the pulse signal from t„ to t^). The weighting 
factor may be defined by the time lapse between the two events n and n+L One possible 
disadvantage of this method is that the system needs good timing resolution because the 
weighting factor is sensitive to time resolution — especially with a small time lapse. 
Additionally, it may need a real-time calculation or a look-up-table to generate the result. 

15 Weighted digital integration 

Digital integration means the instantaneous signal is continuously sampled, and the 
samples of each event are accumulated in digital form. Unlike analog integration, there is 
no dead-time of integration discharging for digital integration. However, this method 
requires a fast analog to digital convertor and a digital accumulator. The rest of the 
20 calculation for estimated total energy is the same as described above — the accumulated 
value may be weighted by a factor that is defined by a time lapse of two countinuous 
events. 

Remnant correction 

Remnant correction is an important step in creating a pileup-free signal, especially 
25 for multiple/continuous pileups. In remnant correction, the residual signal is subtracted 
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from the total estimated energy. The residual term may not be a small correction, and 
without the residual-subtraction, the weighted value may be larger than the real energy of 
the present event by two times or more since the preceding pulse can itself be riding on 
several orther earlier signals, which makes the proceeding weighted value much larger than 
5 the energy of the preceding signal and the energy of the present signal. 

FIG. 22 illustrates one embodiment of remnant correction, where T n+1 is the time 
lapse between event (n+2) and event (n+1), (T n+1 - t n+2 - 1^). A remnant signal may be 
created by weighting the S n+1 with a factor exp(-T n+1 / x\ which may be given by a look-up- 
table (LUT) of T n+1 . This remnant signal may be stored in a register for the residual 
1 0 correction of the next event. 

If the total energy (weighted value) S n is generated from the integrated value I„, 

S = / where I = \ -*-e Ht ~ ttt)/T dt 

From equation (2.1), we have 

1 e~ TJr 1 1 

En+l = l-e~ T ** fT In+l ~\-e~ T * (T In = l-e~ T ^ /T Ift+] ~ ^T^ 7 " Where 
15 T n =U-t n . 

FIG. 23 illustrates this. The look-up-table (LUT) gives the weighting factor of 

1 

t ^ TJr for the input of T n . All the remnant signals (from multiple pileups) that goes to 

the next event may be generated by subtracting the integrated value from the total estimated 
value. This remnant signal may be stored in a register and used for the residual correction 
20 of the next event. 

For signals with two or more decay time 

Some detectors emit signals that follow two or more exponential decay-time 
constants. In these instances, the extractions of pileup-free energy and position may be 
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more complicated. The same general methodology of the common mono-exponential case 
as described above, however, may be used. Each integration may be sampled twice (or 
more) over the cause of integration. The additional sampling provides another set of 
equations to derive the relative magnitutes of the two (or more) exponentials in the total 
5 energy (present event plus remnants). Even though for an isolated event (non-pileup), the 
relative magnitutes for the event are constant and known, the remnant energy in a pileup 
case will have different relative magnitutes for the two exponentials depending on the time- 
lapse between consecutive events. The additional samplings allow one of ordinary skill in 
the art to derive the relative strengths of the exponential. 

10 High resolution time measurement 

High resolution time measurement may be an important issue for pileup processing. 
One needs to measure the time lapse between every two continuous events for the remnant 
corrections of both energy and position signals: see equations (2.1), (2.2) and (2.3). And, 
one also needs the time lapse information to calculate a total energy signal inside the 
1 5 scintillation if the the total value is calculated from the integrated value by equation (2.5). 

Suppose the time laspe measurement error is 8. Then, the measured remnant 
correction value will be is S n * e ~ {t ^ +S)l7 = e~ 5!x * (S n * e'^ 1 ^) . It is e 5/T times of the 
real remnant value. If 5 = 10ns and x = 200ns, then one has e 6/T = e " 10/200 =0.95. That 
means the measured remnant value is 5% less than the real remnant value, and this 5% 
20 remnant vaule still remains in the corrected output result. In a high count-rate situation, 
since the average time lapse between two events is small and the remnant value is big, the 
5% error of remnant value becomes significant. 

This disclosure involves methodology to measure the time lapse with a resolution, 
in one embodiment, less than 1 ns. This means that the remnant correction error may be 
25 less than 0.5% for t = 200ns. A counter with a regular resolution (for example 12ns) 

combined with a high resolution delay line (for example, Ins resolution) may be used to 
measure the time lapse in high resolution. One could use a Ins clock to count the time 
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directly, but it would need faster ECL circuits, which are expensive and may increase digital 
noise level. A better method is described below. 

An embodiment of this disclosure is to use a slow clock to measure the time lapse to 
get a raw time T m . In FIG. 28, T m is 7 cycles of the slow clock. The rising edge of each 
5 event changes the status of a 'flag' signal (from 'high' to 'low\ or from 'low' to 'high'). 
This flag signal inputs into a delay line with multiple tap outputs, which has, for example, a 
1 ns time delay difference between two continuous taps. The rising edge of the first clock 
after each event-trigger latchs the delay line outputs. To measure the time Tc^ between the 
rising edge of the 'flag' signal and the first clock signal, one may count the number of delay 
10 line taps with 'high' output. This number multiplied by the tap delay step (Ins) is the time 
Td n (10nsinFIG. 28). 

The same methodolgoy may be used to measure the time Td I1+1 between the falling 
edge of the 'flag' signal and the following first clock — one may count the number of delay 
line taps with 'low' output and multiply by the tap delay step. In FIG. 28, this is 3ns. The 
1 5 time lapse T, or (t^- tj, between event (n) and event (n+1) may be calculated by T = 
T^Tdjj -Td n+1 In this way, the measured T has the same timing resolution as the delay 
line's, which is a Ins resolution. 

In general, this methodology involves determining an approximate time between the 
two consecutive events using a first clock. Then, a first and second correction to the 
20 approximate time is determined using a second clock, the second clock having a finer time 
resolution than the first clock. Finally, the time between two consecutive events is found by 
applying the first and second corrections to the approximate time. 

With the benefit of this disclosure, those having skill in the art will recognize that 
this methodology will apply to clocks having different clock resolutions {i.e., other than 1 
25 ns example given above). 
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Ultra High Count-rate Application 



To further increase the count-rate in a large detector with many photomultiplier 
tubes (PMTs), multiple HYPER zone circuits may be implemented, A detector may be 
divided into a plurality of smaller HYPER zones by grouping PMTs into multiple smaller 
5 Anger cameras in electronics. Some PMTs may be shared by two or more HYPER zones. 
Events collected by shared PMTs may trigger two or more HYPER circuits at the same 
time. 

This disclosure provides a method to identify a zone that is hit by a real gamma and 
reject false hits within other zones. FIG. 25 shows a detector divided into four zones (zones 
10 K, L, M and N). Since every two neighbor zones share one row or one column of PMTs, 

one gamma hit may trigger two zones if its energy is collected by the shared PMTs. To find 
out the zone that has the "real" hit, a threshold may be set for the pre-normalized position 
signals of each zone to reject the false hit from its neighbor zone. 

FIG. 26 shows an example of a two-zone design. Each zone has one independent 
15 block detector coupled to 4 PMTs. The pre-normalized position of zonelis X1=A+D, and 
the pre-normalized position of zone2 is X2=C+F. If a gamma hits zone2, there is no 
contribution of light to XL However, if a gamma hits zonel (even in the corner crystal 
located in the center of PMT-B or PMT-E), it may still have, in one embodiment, more than 
about 3-5% of total energy going to PMT-A and PMT-D. That means that the pre- 
20 normalized position signal (A+D) will have more than 3-6% of total energy of a real hit. 
Therefore, in data acquisition, a "real" hit in a zone should pass both the energy-acceptance 
threshold and pre-normalized position-acceptance threshold, which in this example would 
be about 3-6% of E0 (the zonal energy-acceptance threshold). However, a false hit caused 
by its neighbor's zone may be accepted by the energy threshold but it will not be accepted 
25 by the pre-normalized position threshold that is set to 3-6% of E0. 

The simulataneous acceptance of the energy threshold and the pre-normalized 
position threshold is one way to reject the false trigger (event) within a zone. This 
technique leads to very robust discrimination because it is immune to differences in 
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electronic-gain balance between different zones. This design also reduces electronic 
complexity and cost. 

Coincidence Application 

The HYPER method combined with multiple zones may significantly increase the 
5 count-rate capability and dose efficiency of PET and gamma cameras. But for coincidence 
detection applications and other timing-detection applications such as positron coincidence 
imaging, a special additional technique is provided. In conventional methods, the signal 
integration time for each event is fixed. Therefore, the end of integration and signal 
digitization (for energy and position) is synchronized with the trigger signal (delayed by the 
10 fixed integration time). The trigger signal is generated at the leading edge or the arrival 

time of the event. Coincidence detection may be performed between the end of digitization 
signals or the trigger signals. 

HYPER, however, uses a dynamic integrating method, which means that the 
measured energy and position signals are no longer synchronized with the leading-edge of 
1 5 the present triggering signal because the energy and position signals are generated or 

digitized at the arrival of the next event (arriving at a random time from present event). To 
detect a coincidence event using the HYPER method, the original trigger signal may be 
delayed by a fixed time, and a new synchronization process may be setup between the 
delayed trigger and the energy/position signals before they go to a coincidence procedure. 

20 In one embodiment, the synchronization process works as follows. The energy and 

position signals of all events may first be buffered into a FIFO (first-in-first-out) temporary 
memory. Then, a delayed trigger signal may be used to read them out individually. Since 
the readout of the energy and position signals from the FIFO queue is time-latched by the 
delayed trigger, and since the delayed trigger signal is synchronized with the trigger signal 

25 (except for the fixed delay-time shift), the energy and position signals may be synchronized 
to the trigger signal. This is shown in FIG. 27. 
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The fixed delay time should be longer than the maximum integration time. The 
trigger-delay circuit should have a good time resolution and a small dead time since the 
delayed trigger signal is synchronized with the trigger signal (except for the fixed delay- 
time shift). 



5 The following table gives an example of processing five continous events by using a 

FIFO and the delay technique of this disclosure to restore the synchronization between the 
event trigger signal and the data output signal (except for the fixed delay-time shift). 





Event number 


1 


2 


3 


4 


5 


Q 


time of event arrival 


0 


0.2 


0.7 


1.5 


2.0 


m 


event integration 
(sampling) time = input 
time of event into FIFO 


0.2 


0.7 


1.5 


2.0 


3.0 


0 


Event readout time from 
FIFO 


1.5 


1.7 


2.2 


3.0 


3.5 




Waiting (holding) time of 
event in FIFO 


1.3 


1.0 


0.7 


1.0 


0.5 



In the above table, the units for time are microseconds. The event integration 
10 (sampling) time is the next event's arrival time if the current event is piled-up by the next 
event. If there is no pileup on the current event, it will be integrated for a fixed maximum 
time period. In this example, 1 microsecond was used. The waiting time is the time 
difference between the FIFO readout time and the event input time. 

The order of the data in a FIFO is the same as the order of events arriving. Hence, 
15 the linear structure of the FIFO is used to line up all the events in the order that they 

arrive. Then, a large fixed delay (larger than the longest integration time) is applied for 
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the event-trigger, and the delayed event trigger is used to "clock" the readout of events 
from the FIFO. Hence, each event in the FIFO may be clocked out (readout) 
synchronously with the event triggering time (exact for a fixed delay). It is the waiting 
(holding) time of an individual event in the FIFO (as controlled by the delayed trigger) 
5 that synchronizes the FIFO readout time to the original event-arriving time. 

The following examples are included to demonstrate preferred embodiments of 
the invention. It should be appreciated by those of skill in the art that the techniques 
disclosed in the examples which follow represent techniques discovered by the inventor 
to function well in the practice of the invention, and thus can be considered to constitute 
10 preferred modes for its practice. However, those of skill in the art should, in light of the 
present disclosure, appreciate that many changes can be made in the specific 
embodiments which are disclosed and still obtain a like or similar result without 
departing from the spirit and scope of the invention. 

WORKING EXAMPLES 

15 In studies with an exemplary embodiment, the present invention was used in 

connection with a 3"x 4" Nal(Tl) scintillator with a photomultiplier from Ortec 
Corporation (Oak Ridge, TN). This detector was connected to (1) regular detector 
electronics setup using a pulse-shaping amplifier (0.5 jas shaping time) and a 
multi-channel analyzer, and (2) to the pile-up prevention circuit ("PPC") of the present 

20 invention. 

The circuit may be interfaced to a computer. For example, the present invention 
was connected to a PC computer with a 133 MHz PENTIUM processor and a high-speed 
input/output board for the studies. The data acquisition software was written with 
LAB VIEW (National Instruments, Austin, Texas). A 99m Tc point source (140 KeV 
25 gamma ray) in air was used for all the count-rate studies. The pile-up-prevention circuit 
of the present invention may also be used to measure the number of nonpile-ups, single 
pile-ups, and multiple pile-ups. 
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For comparison studies, circuits implementing three methods were setup and 
tested. The methods included: (i) the pulse-shaping method (500ns); (ii) delay-line 
clipping (256 ns); and (iii) the method according to the present invention. The 
pulse-height spectra, energy resolution and acquired true-count fractions were collected 
5 as a function of true-even rates using a 99m Tc (140 KeV) point source, which decays with 
a 6 h half-life. This half-life was used to calculate the true-even rate impinging the 
detector when very high count-rate data were acquired. 

EXAMPLE I 

The pulse height spectra of 99m Tc (140 KeV gamma-ray) for the pile-up prevention 

10 circuit of the present invention and the pulse-shaping method (500ns) are shown in 

FIGS. 12A and 12B. FIG. 12A show results of the pulse-shaping method, and FIG. 12B 
show results of the present invention. At 50 Kcps, both measured 10.9% energy 
resolution. At 120 Kcps, the pulse-shaping method (FIG. 12A) started to demonstrate 
pile-up at the higher energy side and distortion at the lower energy side. Above 

15 200-300 Kcps, the spectra were not usable. The PPC method spectra (FIG. 12B) 

maintained the spectra shape even at 2 Mcps. This first study demonstrated that the 
method of the present invention can significantly extend scintillation counting rates. 
However, the energy resolution was poor at 2 Mcps for this first study. Part of the reason 
was due to DC-level instabilities because the resistor divider of the photomultiplier 

20 (PMT) in the ORTEC probe was AC-coupled and the biasing-current was not designed 
for such a high count rate (there was no reason to design for such high rates previously 
because of Nal(Tl) pulse pile-up). Another reason was that the prototype board 
incorporating the apparatus of the present invention was built on a simple board with 
holes and the wiring was done with wirewrapping wire. A third reason was that a single 

25 integrator was used, which may not have been sufficiently discharged if the last event had 
a very high energy. Both the PMT divider and the prototype circuit were subsequently 
improved. Dual integrators and ADC channels ("ping-pong") were provided to assure 
complete discharge of the integrators. The pulse height spectra were measured again for 
the method of the present invention, and the energy resolution was significantly 
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improved, as shown in FIGS. 13A-D. Even at 2 Mcps (FIG. 13C), the energy resolution 
achieved 15% compared to 10.9% at 50 Kcps; part of the reason for this difference is the 
shorter average signal collection time when the count rate is very high. At lower count 
rates, the integration time for most pulses were integrated to ljas, while at 2 Mcps, the 
5 average time-lapse between events was 500 ns which implies that the average signal 
collection time would be equal to or less than 500 ns. 

EXAMPLE II 

In a second study, the method of the present invention was compared with a 
delay-line clipping (DLC) experimental setup (using NIM electronic modules and a 

10 256 ns clipping). The energy resolution results are shown in FIG. 14. At lower count 
rates, the PPC method had better resolution than the DLC method, while at 2 Mcps the 
energy resolution was about the same. The study also generated the photopeak-count 
rates, .i.e. the non-pile-up good counts acquired by each method using a narrow energy 
window, as a function of true-count rates, as shown in FIG. 15. For the method of the 

15 present invention at 2 Mcps, the photopeak counts were 85% of the total incident counts, 
while the DLC had only 25% good counts. The low photopeak fraction in the DLC 
method was due to a high 75% count-loss due to pile-ups which placed the pile-up-event 
energy above the photopeak acceptance-window. The DLC method delayed the 
occurrence of pile-up due to its shorter pulse width, but at very high count rates, its fixed 

20 pulse-width still made it more vulnerable to pile-ups (FIG. 16) than the PPC method, 
which has a variable signal-integrating time and count-recovery capability. Hence, the 
photopeak counts collected by the PPC method is 3.3 times more than the DLC method at 
2 Mcps. The pulse-shaping method only maintained 35% at 0.25 Mcps, and there was no 
useful photopeak data above 0.25 Mcps because of pile-ups (FIG. 15). 

25 The importance of an effective remnant light subtraction for every triggering 

event is demonstrated in FIGS. 17A-F, especially at very high counting rates. At 
600 Kcps, a simple weighted-sum method (SWS) without remnant subtraction yielded an 
acceptable but worse spectrum compared to the present invention. However, at high 
counting rates 1-3.5 Mcps, the energy measurement errors of SWS are significant. Above 
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2 Mcps for Nal(Tl), most events are riding on the signal of other events (see Table 2 and 
FIG. 16), and SWS would measure large energy errors on most events and cause severe 
count-loss outside the photopeak acceptance window. The SWS energy measurement 
errors demonstrated here (shaded spectra in FIGS. 17A-F) should be similar to that shown 
5 in U.S. Patent No. 5,430,406 discussed above. FIGS. 15-17 demonstrates the count and 
energy recovery/restoration capabilities of the present invention. 

EXAMPLE III 

As shown in FIGS. 18A-C, a Monte Carlo simulation of the torso (with liver, 
kidneys, and bladder) was performed to demonstrate the effect of pile-ups in imaging at 

10 2 Mcps for: (a) the method of the present invention (FIG. 18 A); (b) the pulse-clipping 
method (256 ns) (FIG. 18B); and (c) the regular fixed integration method (1 (as) (FIG. 
18C). The study simulates the electronic processing of each method in a gamma camera. 
The results are shown in FIGS. 18A-C. At 2 Mcps, the present invention method showed 
all the uptake organs; the conventional method shows no organs; the pulse-clipping result 

15 was in between the two, with high pile-up artifacts between the two kidneys and 
diminishing activities in the bladder and kidneys. 

Besides preventing measurement errors from pile-ups, the present invention 
allows data acquisition time to be shortened significantly. This fast data acquisition is 
important in commercial/medical applications, as data acquisition time can be 

20 significantly reduced. As discussed above, the count loss with the conventional method 
is very high above 160,000 cps and saturates at 200,000 cps. Although the data appears 
to indicate that the conventional method can approach 200,000 cps, the pile-up fraction at 
this limit is very large, and the photopeak is badly distorted. Therefore, a large 
percentage of the data is not useful The present invention can count over 2,000,000 cps 

25 with very little count loss, while preserving a very good photopeak, and all of the data is 
useful. Thus, the present invention is a significant improvement over the prior art. 

Data from an embodiment of the present invention demonstrate that the invention 
significantly improves measurement of the energy of gamma rays at very high 



25075668.1 



-53- 



count-rates, even when multiple pile-up approaches 60% and single pile-up approaches 
17% (nonpile-up at only 13%) at 2,000,000 cps. The data indicate that a 10-20 fold 
improvement in count-rate capability is feasible with the present invention. 

These studies discussed above showed that the methods of the present invention 
5 are feasible. Compared to the pulse-shaping method, it has a 15 times higher counting 
rates, and yet at low counting rates, the PPC energy resolution is the same as the standard 
pulse-shaping method. The PPC energy resolution only degraded slowly from 10.9% at 
50 Kcps to 15% at 2000 Kcps. Compared to the pulse-clipping technique, the PPC 
method has better energy resolution at regular count rates (10.9% vs. 15%) and 3.3 times 
10 higher photopeak (non-pile-up) counts at 2000 Kcps. The count and energy recovery 

capability of the present invention at very high count rates is important for extending the 
maximum count rates of scintillation detectors. 

All of the methods and apparatus disclosed and claimed herein can be made and 
executed without undue experimentation in light of the present disclosure. While the 

15 compositions and methods of this invention have been described in terms of preferred 

embodiments, it will be apparent to those of skill in the art that variations may be applied 
to the methods and apparatus and in the steps or in the sequence of steps of the method 
described herein without departing from the concept, spirit and scope of the invention. 
All such similar substitutes and modifications apparent to those skilled in the art are 

20 deemed to be within the spirit, scope and concept of the invention as defined by the 
appended claims. 

APPENDIX 

(1) Energy Pile-up Prevention Formulae Algebra 

Suppose at time t j9 a y-th gamma ray (y } ) is detected. The j-th gamma ray is a pile- 
25 up event on top of two preceding gamma rays y yi and y^. From Knoll, 1979, the 
instantaneous scintillation emission of a scintillator after detecting a gamma ray is 

m-'e**'" (i) 
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where E is the total scintillation signal generated by the scintillator, and E is proportional 
to the energy of the gamma ray detected. The instantaneous signal at time tj for 
tj < t < tj +I , contains emission from three gamma rays y y2 , y^, y j? and is given by, 



q (t) = iL e <»,W + !j± e <"»>h +El± e ^,)h 

XX X 



E 



n <t-t oh 



^h 



(2) 



I «j-'jJh 



(3) 



If q(t) is amplified or weighted by an amplification factor t: 



(4) 



10 



If q(t) is only integrated from t, to t (this implicitly means that the integrator was 
discharged or zeroed just before t,) from equation (3), 



| g (u)du = tesE^'^+E^'^-e**] 



(5) 



If equation (4) and (5) are summed, the weighted sum within the time period of gamma- 
ray yj will be: 



Sj s *ft)+ EjAej-^Ej^^ 



-{'j-tjJh 



(6) 



The same analysis is now applied to the preceding time period where < t < tj (the time 
15 period of y H ). From equation (1) and FIG. 13, the instantaneous signal after the onset of 
Ym is: 



(7) 
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The weighted sum for this period is: 



Hence is the term inside the square bracket of equation (6), and equation (6) 
becomes: 



5 SrEj+Sj-^' 1 " 1 ^ 



or 



Ej = Sj-S M e [tj ~ >J - >h 



(9) 



Hence the energy deposition of the y-th gamma ray can be derived from the y-th 
weighted sum minus the preceding weighted sum decreased by a emission decay factor, 
provided that the integrator is discharged immediately before the arrival of the y-th 
10 gamma. Note that the weighted-sum Sy as shown in equation (6) is not a measurement of 
energy of the y-th event, it is the measurement of the total energy trapped inside the 
scintillator just after the y-th event hit; i.e., S } includes the energy of the y-th event and all 
the remnant energy of previous multiple-pile-ups when the y-th event is deposited. 

Equation (9) is an exact solution, regardless of how many previous gamma tails 
15 on which the y-th gamma is riding. The present invention thus provides the first method 
that attempts to recover all triggering events and compensating for multiple-pile-up 
exactly. Thus, the present invention may achieve higher count rates than possible with 
known methods. 

Equation (9) can be easily understood intuitively. If there is no pile-up, Sj is the 
20 detected energy of the y-th gamma ray (independent to when Sj is sampled). If there are 
pile-ups, Sj is the total energy in the detector in the y-th period (independent to when Sj is 
sampled), which includes the detected energy of the y-th gamma ray and the remnant 
emission of all the preceding gamma ray. This y-th remnant emission is simply the total 
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energy in the preceding period (S H ) decreased by the emission decay factor governed by 
the time difference between the y-th gamma and its preceding gamma. 



(2) Pile-up Prevention Formulae for Gamma-Camera Anger-Positioning 

In Anger-position method, assuming that no pile-up occurs, the y-location is 
5 estimated by the centroid calculation: 



where Ef is the total signal received by Photomultiplier-i (PMT-i) when a y-ray is 
detected, and is the physical location of PMT-i. X is the prenormalized pseudo- 
position signal which needs to be normalized by the total detected energy Z of the y-ray, 
10 to generate the centroid location X. X can be defined as; 

X' = X + -X„ 

If at time t\, the first y-ray with energy Z; is detected, the total energy (after signal 
integration) distributed to PMT-1, PMT-2, . . . , PMT-n are E] j,E]2,..., Ej n , 
respectively. The total integrated energy detected by the camera is: 



(10) 



15 



li 



(ll) 



From equation (10), the prenormalized position X for conventional Anger 



positioning is: 



(12) 
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In conventional Anger-position, En in equation (12) are the integrated PMT 
signals (integrated for 0.6-1 us), and therefore^ in equation (12) is a static quantity 
obtained after a fixed 1 jus integration time. Hence, if there is a pile-up during the 
integration time, X x will be wrong. In the method of the present invention, the 
instantaneous prenormalized position X(t) is used: 

*.'(') =E i ^e < '-" ,A 

The weighted instantaneous prenormalized position is: 

tf;(')=Z,*Ae ,, -" )/ ' (B) 

The integrated prenormalized position is: 

Ix^u^^uk-e^^) (14) 
If equation (13) and (14) are summed, the weighted Sum of X j(t) is given by: 
SXl = xX[(t) + { X[{u)du = ^X t E u (15) 

From equation (12) and equation (15), we obtain SX[ = X[ . Hence , if the event 
yj is not a pile-up on previous events, the weighted sum, SX[ , is also an estimation of the 
Anger centroid-location for the gamma ray (yj, just like static quantity X[ in equation 
(12) for conventional Anger-method. Furthermore, the weighted sum SX[ is a time 
independent quantity which can be sampled at anytime and still provide an unbiased 
estimation of the prenormalized location. Therefore sum position SX[ is just like its 
energy-counterpart weighted sum Sj in equation (6). 

However, if the event (ft) is a pile-up on a previous gamma-ray (y 0 ), the 
prenormalized position weighted sum will not be equal to the true prenormalized position 
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of y l5 because SX[ contains the remnant prenormalized position weighted-sum SX' 0 
generated by the remnant scintillation of y 0 . The y 0 remnant position weighted-sum (a 
position vector) adds an error-position vector to the y ] position. The magnitude of this 
error vector increases if the time lapse (t r t 0 ) between y 0 and y x decreases. These 
5 positioning errors are especially severe when the count rates are very high that (a) most of 
the events ride on the remnant of previous events and (b) multiple-event pile-ups are the 
norm. Hence, the general solution to derive the true prenormalized position of each event 
Yi requires the subtraction of the remnant-position signals of previous events from the 
prenormalized position weighted-sum, SX[ and SYj , triggered by each y { . 

10 In the general case where any event can be riding on the signals of one or more 

previous events, the weighted instantaneous prenormalized position signal triggered by an 
incident event y m is given by the superposition, 

*K^Y.x^e M (i6) 

where t m < t <t m+1 . The corresponding integrated prenormalized position signal is given 
15 by: 

^Z^-ue^H-e^^j (17) 



This integration can be carried out by first zeroing the integrator just before the 
onset of Y m , and then restarting the integration immediately. 

Summing equations (16) and (17), the weighted sum of the prenormalized 
20 position signal is: 
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i b (18) 



It can be observed that the term inside the bracket in equation (19) is simply SXm - 1 (see 
equation (18)). The recurring relationship reduced equations (19) and (12) to 

Therefore, the prenormalized Anger-logic position (Xm , Ym ) of y m can be derived from 
the weighted prenormalized position sums of y m , y m _ 1? and the time-lapse between the two 
events, 

X'm= SX'm SXm - 1 g'^-'^ 

10 Since the derivation of equation (20) assumes that y m is riding on any numbers of 

preceding events, equation (20) provides an exact determination of the position of y m , 
even when it is part of a multiple-event pile-up. 

The pile-up-free prenormalized position ( X'm, Y'm ) will then be renormalized by the 
energy of y m as in conventional Anger-logic to remove the energy- scaling effect in 
15 generating the position of y m . Hence, the energy of y m also has to be extracted correctly 
from a multiple-event pile-up if the position of y m is to be decoded accurately. The pile- 
up-free energy of y m can also be derived from the energy weighted-sum of y m , y m .j, and the 
time lapse between these two events, 
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E m =S m -S m _, e"^> (21) 

where S m is the sum of the weighted instantaneous scintillation signal and the integral of 
the instantaneous scintillation signal. 
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